What determines the spreading of a wave packet? 
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The multifractal dimensions D% and of the energy spectrum and eigenfunctions, resp., are shown 
to determine the asymptotic scaling of the width of a spreading wave packet. For systems where the 
shape of the wave packet is preserved the fc-th moment increases as t hli with /3 = D^jD^, while in 
general t kf3 is an optimal lower bound. Furthermore, we show that in d dimensions asymptotically 
in time the center of any wave packet decreases spatially as a power law with exponent — d and 
present numerical support for these results. 
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For the case of a free particle, the spreading of a quan- 
tum mechanical wave packet is a textbook example. One 
finds that asymptotically the width increases linearly in 
time and the corresponding energy spectrum is absolutely 
continuous. In the case of a point spectrum, on the 
other hand, there is asymptotically no increase of the 
width of a wave packet, if the eigenfunctions are semiu- 
niformly localized Q. Between these extremes flourishes 
the world of quantum systems with fractal energy spec- 
tra and eigenfunctions. They include systems studied in 
the early days of quantum mechanics, like Bloch elec- 
trons in a magnetic field || as well as quasicrystals jjj 
and disordered systems, e.g. at the Anderson transition 
H or in the Quantum Hall regime || . A natural question 
then arises: What determines their long-time dynamical 
properties, e.g. the spreading of wave packets (quantum 
diffusion) and the decay of temporal correlations? 

Very few rigorous and general answers are known so 
far. Temporal correlations decay as t~ D % [||, where 
is the correlation dimension of the spectral measure /j, 
(i.e. the local density of states). For the spreading of 
a wave packet, on the other hand, two inequalities were 
derived by Guarneri. The growth of the fc-th moment 
m k (t) ~ t kf3k is bounded from below by D? < (3 k @ and 
the entropic width grows faster than i D i || , where is 
another generalized dimension of the spectral measure. 

In search for equalities, relating the growth of the mo- 
ments of a wave packet to fractal properties of the spec- 
trum, numerical studies of the Harper model P,|i"o||, the 
Fibonacci model (l^Jl^], and the kicked Harper model 
jl3| were performed, suggesting, e.g. $2 = D , in agree- 
ment with a heuristic argument f^Tcl, where Dq is the 
fractal (box-counting) dimension of the spectrum. More 
recently, numerical studies |l4|,[l5| showed that this sim- 
ple relation does not hold exactly. These studies also re- 
vealed that one often finds multiscaling in time |lq] , i.e. 
Pk varies with fc. For a restricted class of systems, Man- 
tica showed that (5k = D^_ k , whereas in general it is at 
best approximate pffl. The similar relation (5k — -Di-fc, 



now with the multifractal dimension of the ("global") 
density of states, was proposed to hold after averaging 
the dynamics of all wave packets started at different ini- 
tial sites jUI. Therefore, the basic question, of what 
determines the spreading of a single wave packet, still 
remains open. 

In this paper we derive a partial answer to this question 
(Fig. 1), by including, for the first time, fractal properties 
of the eigenfunctions. For systems where asymptotically 
the fc-th moment of a wave packet increases proportional 
to t k $ , we show 
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where D% is the correlation dimension of the (suitably 
averaged) eigenfunctions. For the more general case of 
multiscaling in time we derive an optimal lower bound 
for positive moments under reasonable assumptions: 



(3 k > D$/D. 
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FIG. 1. Quantum diffusion of a wave packet initially started at 
site no- The staying probability decays as t i, whereas the fc-th 
moment increases as t fe ft . If the form of the wave packet remains 
unchanged (uniform scaling) Eq. (|lj) leads to (3^ = /3 = Dj/O^ 
for k > 0. For the case of multiscaling dynamics we find the lower 
bound f) < for fc > (Eq. (Q)). In any case, the center of the 
wave packet spatially decreases as \n — no\ D 2 _1 for | n — ng | <C ip . 
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We verify these results numerically for the on-site Fi- 
bonacci model and the Harper model, where it is a much 
better lower bound than the previously obtained D± , 
which did not make use of fractal properties of the eigen- 
f unctions. Surprisingly, we find from these dynamical 
properties that in d dimensions asymptotically in time 
the center of any wave packet decreases spatially like a 
power law with an exponent — d (see Fig. 1). 

The moments of a wave packet tp{n,i) 1 initially lo- 
cated on site no, are given by rrik{t) — J2 n ^n \ n ~ 
no\ k \(p(n, t)\ 2 , where for convenience we use a finite one- 
dimensional lattice. To derive the scaling properties of 
the moments we need a relation between the integrated 
wave packet 



amplitude. Assuming a power-law scaling one there- 
fore has a function G of a single variable defined by 
G(x/tP) = g{x,t). This is equivalent to the property 
that the positive (k > 0) moments scale as m,k{t) ~ £ fc/3 , 
which we call uniform scaling. From Eq. (||) and the 
properties of G we find that S(x, ui) is a function of xuj 13 
only. Together with the scaling behaviors of Eqs. M) 



and 



one easily finds the result (3 — / D\ 



(Eq. 



From Eq. (|^) we can thus even determine all negative mo- 
ments: nik(t) ~ t kf3 for k > — n ' 



-D%, whereas mjfe(i) ~ t 



i(x,t)= Yl MMI S 



\n—no | <x 



and the spectral function [fjj 



for k < -D$. 

Secondly, we now want to study the case of multiscal- 
ing dynamics. We cannot do this in the most general 
form, but need two reasonable assumptions. We first 
assume that for xu>P <C 1 the spectral function S(x, to) 
(3) is determined by its limiting scaling behaviors (Eqs. (Q) 
and (§)), namely S(x,u) ~ x D * ui D * §§. Using Eqs. (§) 
and (6) this implies uniform scaling for the center of the 
wave packet only, i.e. g(x,t) = G(x/t 13 ) for x <C t" . In 



S(x,w) = V V (no)^i(«o)^i(n)^J(n) fact ' {t follows 9fat) ~ x ° 2 1 ° 2 for x < ^ ■ Remark 
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ably, this even determines the shape of the center of the 
wave packet, namely, an algebraic decay |n — tiq\ D 2 _1 for 
\n — no | <C ir . In addition, we assume that outside this 



with Ei and tpi denoting eigenenergies and eigenfunc- 
tions, respectively. For systems described by a real sym- 
metric matrix these are simply related by 



region it decays faster than 
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D., 



which it cer- 



and 



g(x,t) — / duicosuit— — S(x, ui) 
du> 



, 2 f°° , sinwf , . 
b(x,ui) = - dt——g{x,t). 
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For the special case x = 1 the spectral function simpli- 
fies and we assume the following scaling behavior pQ|: 



S(x = 1, ui) 



\Ei-Ej\<u 



and 



(8) 



where ui m [„ is the smallest energy scale where the scal- 
ing of Eq. (0) still holds for a given finite system. Physi- 
cally, u>rmn 1S ^ ne mverse °f the time when a wave packet 
reaches the boundary and when the power-law moment 
growth is modified. Drj; is the (standard) correlation di- 
mension of the spectral measure and D2 denotes the cor- 
relation dimension of the averaged eigenfunctions j|21f| . 

We will first study the special case that the over- 
all shape of the wave packet stays the same during its 
time evolution, with, of course, a scaling in width and 



tainly has to do asymptotically in order to stay normal- 
ized. Under these assumptions, it follows that for a large 
enough time T there exists an X such that g(x, T) oc x D i 
for x < X and that for all t > T with g(£(t), t) = g(X, T) 
the relation £{t)/t fj > XlT fi must hold. Using a result by 
Guarneri and Mantica 0, stating that all (3% (for k > 0) 
are larger than the scaling exponent of £(t), it then fol- 
lows that all positive moments scale as t kfik with /3* > [3. 
Therefore (3 is a lower bound on the scaling of the posi- 
tive moments that incorporates multifractal properties of 
the eigenfunctions. In fact, for negative moments (k < 0) 
we find that [3 > (3k is an upper bound. In that sense (3 
is an optimal lower bound for all positive moments. 

Furthermore, note that this analysis can be readily ex- 
tended to systems in higher dimensions d by appropri- 
ately generalizing the definitions of g and S. Under the 
condition < 1 Eqs. ([!]) and (J2j> remain unchanged and 
the spatial decay of the center of the wave packet has an 
exponent — d. Finally, this formalism can easily de- 
scribe the averaged dynamics of wave packets started at 
different initial sites no. To this end we introduce an av- 
erage m Eqs. (|) and (|) defi ning exponents and , 
which replace D% and in the above results pj| . 

One can try to understand the derived results in an in- 
tuitive way: The temporal decay of the center of the wave 
packet is known to be given by t-°2 0. Naively, nor- 
malization then requires a spreading described by an ex- 
ponent [3 — £>2 /d in d dimensions. If the spreading, how- 
ever, takes place in a space with an effectively reduced 
dimension D% instead of d j§, we have (3 = D%/D$. 
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ilarly, the integrated wave packet g(x,t), which usually 
increases initially as x d , here increases with x D ? instead. 
This leads immediately to a power-law decay \n— no\ D ? ~ d 
for the center of the wave packet in the <i-dimensional 
embedding space. 

In the remainder, we want to give examples and nu- 
merical evidence supporting the above analysis . 

For disordered systems at the metal-insulator transi- 
tion in 2 (with strong magnetic field or symplectic sym- 
metry) and 3 dimensions Eq. is already known. These 
systems are simpler in the sense that the fractal dimen- 
sion Do of the spectrum is 1 and it was shown that 

D% = dD% holds HfOH, Wh6re d = 2 ' 3 iS th<3 Spa " 
tial dimension. On the other hand, (3 = l/d follows from 

the non-fractal value of Dq using heuristic arguments, as 
in Refs. |^,|l^,^2), and was numerically confirmed for the 
second moment ||^7|. The prediction that the asymp- 
totic shape of the wave packet follows the power law 
\n — no\ D 2 ~ d waits to be observed in disordered systems. 

We now consider the on-site Fibonacci chain, which is 
a one-dimensional model of a quasicrystal || given by a 
tight-binding Hamiltonian where the on-site energy V n 
takes the values +V and — V, arranged according to the 
Fibonacci sequence. Figs. 2 and 3 show that the ratio 
(3 = D2/D2 — 0.477 coincides, within the numerical 
accuracy, with the scaling behavior of the moments, thus 
confirming Eq. (Q) g|. 
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FIG. 2. The moments m^{t) for k = 1,2,..., 10 of a wave 
packet initially started at a symmetric site of a Fibonacci chain 
with V = 2 and 17711 sites. The insets show the determination of 
= 0.072 (left) and D% = 0.151 (right) from the limiting scaling 
behavior of S(x,u>) yielding /3 = D^/D^ = 0.477, which is much 
larger than the information dimension = 0.153. 
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FIG. 3. The scaled moments rrife(t) /t k @ of a Fibonacci chain 
corresponding to the parameters of Fig. 2 are (on average) constant, 
thus confirming Eq. (Ill) . 




FIG. 4. The spectral function S(x,u>) for the Harper model 
with a = 6765/10946, an approximant of the golden mean, and 
v = of a wave packet initially started at n = 0. It shows the 
scaling behavior ~i !u J for small x and lo, as we have it as- 
sumed for deriving Eq. (Q). (Deviations close to S(x,lu) = 1 might 
be the cause for multiscaling in time.) We find = 0.134 and 
Df = 0.291. 
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FIG. 5. The values (3^ (diamonds) of the power-law scaling mo- 
ments m k (t) ~ t h $k of the Harper model with cr = 10946/17711 
and v = of a wave packet initially started at n = 0. One can see 
that f3 = / D 2 = 0.460 is a much better lower bound for all pos- 
itive moments than D 1 ^ = 0.227. For negative moments, k < — 1, 
we find k/S^ = —D^- The inset shows the power-law behavior of 

the integrated wave packet g(x,t = 3 X 10") ~ x 2 corresponding 
to a power-law decay \<p(n, t)| 2 ~ \n — no\ D 2 ~ 1 . 
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In general, though, quantum systems with a fractal en- 
ergy spectrum generate multiscaling dynamics p4] ^16[. 
As an example we consider the Harper model ]29| , de- 
scribing an electron in a two-dimensional periodic po- 
tential and a perpendicular magnetic field Qj. For a 
the number of magnetic flux quanta per unit cell it is 
given by a tight-binding Hamiltonian, where now V n — 
2 cos(27r an + v) holds. Fig. 4 shows that the scaling as- 
sumption for S(x, lo) is fulfilled ]3(| and Fig. 5 shows that 
P < 0k is a good lower bound for the positive moments. 
We have verified (3 < fik also for wave packets started at 
other sites no as well as for a few other irrational a. In 
fact, in all cases f3 considerably improves the lower bound 
deduced from spectral properties . In addition, we 
have confirmed the prediction for the shape of the cen- 
ter of the wave packet, namely a power-law decay with 
exponent — 1, i.e. g(x,t) ~ x D 2 (Fig. 5 inset). 
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